Parametric invariant Random Matrix Model and the emergence of multifractality 
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We propose a random matrix modeling for the parametric evolution of eigenstates. The model 
is inspired by a large class of quantized chaotic systems. Its unique feature is having parametric 
invariance while still possessing the non-perturbative breakdown that has been discussed by Wigner 
50 years ago. Of particular interest is the emergence of an additional crossover to multifractality. 
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I. INTRODUCTION 

The analysis of structural changes of eigenstates as a 
parameter x is varied, has sparked aCTeat deal of re- 
search activity for many years 0, S S HB II S S 13 ■ Of 
particular importance are quantized chaotic or complex 
systems where the change of x may represent the effect 
of some externally controlled field (like electric field, gate 
voltage or magnetic flux). Thus, these studies are rele- 
vant for diverse areas of physics ranging from nuclear 
P, |l0| and atomicphysics to quantum chaos ^ |^ 
and mesoscopics 

In all these studies, Random Matrix Theory (RMT) 
played a dominant role as a reference theory that de- 
scribes the universal properties of the eigenstates of com- 
plex systems. RMT was introduced 50 years ago by 
Wigner as a tool to describe the statistical properties of 
the eigenvalues of complex nuclei. Until recently, the ma- 
trices in the frame of RMT were assumed to be homoge- 
neous, i.e., all matrix elements were set to have identical 
statistical properties. Under this simplification random 
matrices are rotationally invariant, a property that sim- 
plifies their theoretical analysis. In physical applications 
this implies that interactions are assumed to be so strong 
and complex that no other parameters, apart from the 
symmetry of the Hamiltonian matrix, are relevant. As 
a result, such random matrices can be associated to the 
extreme case of maximal chaos, which is known to ap- 
pear in various physical systems such as heavy nuclei, 
atoms, metallic clusters, etc. Moreover, one can treat 
full random matrices as a typical model when describ- 
ing local statistical properties of spectra and eigenstates 
in some range of the energy spectrum, typically, in the 
semiclassical region. 

On the other hand, the conventional RMT can not de- 
scribe important phenomena such as localization of eigen- 
states, neither can be directly applied to obtain spectra 
of realistic models. For this reason, much attention has 
been recently paid to the so-called Wigner Band Ran- 
dom Matrix (WBRM) model which is characterized by 
the free parameter 5, that represents the effective band- 
width of a Hamiltonian matrix. Among the important 
applications of the WBRM model we mention the study 



of localization in quasi-one-dimensional disordered sys- 
tems. Also the WBRM model has been applied to the 
analysis of either chaotic or complex conservative quan- 
tum systems that are present in nuclear physics as well 
as in atomic and molecular physics. 

Despite its success, the standard WBRM model has 
severe limitations in modeling realistic systems. We ex- 
plain these limitations in section II, and further motivate 
the introduction of a new RMT ensemble to which we re- 
fer as the Winger Lorentzian Random Matrix (WLRM) 
model. In Sec. Ill, the WLRM model is shown to have 
a "parametric invariance" property which is characteris- 
tic of any realistic system but is missing in the standard 
WBRM model. The analysis of the local density of states 
of the WLRM model is done in Sec. IV. The multifrac- 
tal properties of the eigenstates of the WLRM ensemble 
are analyzed in Sec.V. Our conclusions are summarized 
in Sec. VI. 



II. RMT MODELING 

The pioneering work in this field has been done by 
Wigner , who has motivated the studies of RMT mod- 
els of the type 



E + xB 



(1) 



Both E and B are real symmetric matrices of size NxN. 
The elements of the diagonal matrix E are the ordered 
energies {En}, with mean level spacing A, while B is a 
banded random matrix which is characterized by a band 
profile C(r). Namely, the entries of B are random num- 
bers that are drawn from a normal distribution with zero 
mean and variance given by 



\Bnmn = Cin-m) 



(2) 



For the study of spectral statistics of the energy levels it 
turns out that full matrices (C(r) — 1), say of the Gaus- 
sian Orthogonal ensemble (GOE), are enough in order 
to capture the universality which is found in quantized 
chaotic system. But for the study of parametric evolu- 
tion of the eigenstates it is essential to take into account 
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the band profile, which is dictated by semiclassical con- 
siderations. Namely, C(r) is merely the scaled version of 
a classical power spectrum C(w) which is obtained via a 
Fourier transform of the classical correlation function of 
the generalized force T{t) = —dH/dx. 

A. The standard WBRM model 

The standard WBRM model assumes a rectangular 
band profile: 

For this model Wigner has found that the eigenstates 
undergo a transition from a perturbative Lorentzian-type 
line shape to a non-perturhative semicircle line-shape. It 
should be clear that the Wigner Lorentzian can be re- 
garded as the outcome of perturbation theory to infinite 
order, while the semicircle line-shape is beyond any order 
of perturbation theory. 

The existence of a transition from a perturbative to a 
non-perturbative line shape is a generic feature of any 
realistic (quantized) Hamiltonian. In the latter case the 
semi-circle line-shape is replaced by a semi-classical line 
shape. 

B. The modified WBRM model 

The WBRM model suffers from a serious drawback. 
Unlike generic canonically quantized Hamiltonians, the 
statistical properties of its Hamiltonian matrix are not in- 
variant under x — *■ x -I- const. In fact there exists two wise 
modified versions of the WBRM model HI 13 

which are 

manifestly x invariant by construction. For example we 
cite one of them: 

E + cos(a;)Bi -I- sin(x)B2 (4) 

Here Bi and B2 are uncorrelated banded matrices. It is 
quite easy to be convinced that this model is x invariant. 
One simply has to set x x + const, to expand the sin() 
and the cos(), to define B'l and B2, and to observe that 
they are uncorrelated with the same band profile as Bi 
and B2- 

However there is a "price" for using such modified 
model. It is not difficult to prove that the parametric 
nature of this model is essentially perturbative: The as- 
sociated local density of states does not exhibit the non- 
perturbative crossover that has been highlighted in the 
previous subsection! 

C. The Winger Lorentzian Random Matrix model 

In the present paper we introduce a new RMT ensem- 
ble to which we refer as the Winger Lorentzian Random 



Matrix (WLRM) model. The Hamiltonian is assumed to 
have the standard form of Eq. , and it is characterized 
by the band profile 



There are several good reasons that motivate the intro- 
duction and the study of this model, which we are going 
to clarify: 

(1) There is a major class of quantized chaotic systems 
that can be described using this model. 

(2) Unhkc the standard WBRM model it has the de- 
sired X invariance property that characterizes quan- 
tized models. 

(3) Unlike the common x invariant version of the 
WBRM model it exhibits the transition to a non- 
perturbative line shape. 

(4) The emergence of multifractality, which is absent in 
the WBRM model, is a fascinating issue by itself. 

Let us expand on the first point. We recall that the 
flat band profile of the standard WBRM model is mo- 
tivated by the realization that many observables (say 
J-'{t)) of chaotic systems exhibit "white" power spec- 
trum: Cjrjr{uj) const. However, in many cases it is 
G{t) = T that has the "white" power spectrum 
In the latter case the relation Cjrjr(uj) — Cgg{^)I'^'^ im- 
plies Lorentzian tails. 

III. PARAMETRIC INVARIANCE 

An important feature of a generic canonically quan- 
tized Hamiltonian Ti.{Q, P; x) is its parametric x invari- 
ance. Given x* we can represent the Hamiltonian by a 
matrix 

n = E + 5xB (6) 

where Sx — x ~ x* . If we take two different values of x* 
we get two different B matrices. But if the two values 
of X* belong to the same classically small window, then 
(by definition) the band profile C (r) comes out the same. 
Still from a quantum mechanical point of view a classi- 
cally small range of x values is typically regarded as huge. 
This means that in general, quantum perturbation the- 
ory cannot be used in order to describe the parametric 
evolution within this range. 

In Fig. n we present C(r) for the standard WBRM 
model and for the WLRM model. We see that the profile 
of the perturbation matrix B of the WBRM model is de- 
formed as X is increased, while that of the WLRM model 
remains the same. We have found out that this x invari- 
ance does not hold for other (non-Lorentzian) power law 
profiles. As examples, in Fig. |21 we present the profiles 
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FIG. 1: (a) The band profile C(r) of tlie standard Wigner 
model using b — 500 and x* = 0, 0.1, 1, 3, 7, 15, 31, and 63. 
(b) The band profile C(r) of the Lorentzian Wigner model 
with 6 = 1 and a;* = 2.5, 27, 277, and 2777. The dashed line 
in the inset of (b) with decay ~ is plotted to guide the 
eye. 



FIG. 2: The band profiles for the power law matrix C(r) = 
[1 + (r/b)"]-^ with (a) = 0.5 and (b) 5. 6=1. In (a) 
X* = 2.5, 27, 277, and 2777. In (b) x* = 0.13, 1.4, 14, and 
144. The dashed line in the inset of (a) [(b)] with decay 
~ r"*^'^ [~ r^^] is plotted to guide the eye. 



of B for increasing x in the case of the power-law profile 
C{r) = [1 + [r/bY]^^ with fi = 0.5 and 5. We conclude 
that the x invariance is a unique property of the WLRM 
ensemble. 



IV. LDOS ANALYSIS 

The local density of states (LDOS) is the major tool 
for the characterization of the parametric evolution of 
the eigenstates. The overlap of the eigenstates \n{x)) 
for a given value of x with the eigenstate \m{Q)) of the 
a; = Hamiltonian is P{n\m) — \{n{x)\m{Q))\'^ . This 
can be regarded as a distribution with respect to n. By 
averaging over the reference level m we get the line shape 
P{n — m). Up to trivial scaling this is the LDOS. 

The considerations that are required in order to gener- 
alize the calculation of the LDOS line shape for a general 
band profile have been introduced in ||^|8|. Here, we ap- 
ply such methodology in order to analyze the parametric 
evolution of the LDOS for the WLRM model. 

For X = the LDOS is trivially P{r) — 5^,0 due to 
orthogonality. As x increases, perturbative tails start to 
appear. By employing standard first-order perturbation 
theory we get Pfopt(?') ~ 1 for r = 0, while 



:ir) = 



1 



{En—E„ 



A2 (62 



(7) 



for r 7^ 0. In the second equality we have substituted 



the matrix elements of B using Eq. Q. The above ex- 
pression applies for x < Xc where Xc is the perturbation 
strength needed to mix neighboring levels only. 

For X > Xc Wigner had found, in the frame of the 
WBRM model, that the LDOS line shape can be calcu- 
lated using perturbation theory to infinite order. In case 
of the WBRM model one obtains a Lorentzian. Assum- 
ing the validity of infinite order perturbation theory we 
come out with a Lorentzian-type approximation for the 
LDOS of the WLRM model: 
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A2 (fe2 +r2) [(r/A)2 +r2 



(8) 



Eq. ((SJ is an approximation because all orders of pertur- 
bation are treated within a Markovian-like approach (by 
iterating the first order result) and convergence of the 
expansion is pre-assumed. The energy scale F defines 
the region where a non-perturbative mixing of levels oc- 
curs. Inside this region the perturbative profile P-pwri''') 
does not describe the actual LDOS lineshape. F is de- 
termined by imposing normalization of Ppht('')- For the 
WBRM model it was found l\ that F = a;2C(l)/A. For 
the WLRM model we get 



6A 
~2" 



47ra;2 



(9) 
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FIG. 3: The LDOS lineshape P(r) in the (a) standard pertur- 
bative, (b) extended perturbative, and (c) non-perturbative 
regimes for the WLRM model with 6 = 1, A = 1, and 
N = 5000. For this set of parameters Xc ~ 0.8 and x^^t ~ 1-5. 
The first-order perturbation theory profile Pfopt('") from 
Eq. O is included in (a). The Wigner-type lineshape Pprt(j') 
from Eq. which is expected to be valid in the perturbative 
regime, is also included (no fitting parameters) in (b) and (c). 
The dashed line in the insets of (a) and (b) with decay ^ r"** 
is plotted to guide the eye. The gray line in the inset of (c) 
is a semicircle fitting to P{r). 



Obviously, for F ^ A the (infinite-order) LDOS profile 
^prt('') reduces to the standard first-order perturbation 
theory expression Pfopt('')- Therefore Xc can be deter- 
mined by the condition T(xc) ~ A, leading to 



A 



(10) 



In Fig. 13 we display our numerical results for P{r) 
for the WLRM model with & = 1, A = 1, = 5000, 
and various perturbation strengths x. We see that the 
agreement with the perturbative expression JSJ persist up 
to some perturbation strength Xp^t. Above x^^^ the LDOS 
lineshape P{r) becomes semicircle in complete analogy 
with the WBRM model scenario. 



We want to find the value x^^t up to which the pertur- 
bative expression Pprt(?') describes reasonably good the 
LDOS lineshape. To this end we compare the dispersion 

of Pprt(?-), 



J-BpRT^Ax /^r2PpRT(r) , 



to the dispersion of the actual LDOS 
SE = x^C{r) . 



(11) 



(12) 



Expressions for SEp^j, and 5E in case of our WLRM 
model can be obtained by replacing the sums above by 
integrals: 

SEpp^^x^bib + T/A)-^/^ , 



SE 



2b [7r/2-arctan(l/6)]^/^ 



(13) 



(14) 



Note that for small x, 5£'prt — 5E k x^-KbY/"^. SE is a 
linear function of x for all perturbation strengths while 
for large enough x SEp^^^, becomes sublinear: SEp 



See Fig. ^ The border x^^t is identified as the 



perturbation strength for which dEp^.p{Xp,t) 
where 7 < 1: 



iSEix, 



AVb 



a/tt - 272[7r/2 - arctan(l/6)] 
272[7r/2- arctan(l/fe)] 



(15) 



We typically use 7 — 0.8. 

The vahdity of Eqs. ©, l(T^ . and ((TUl is confirmed in 
Fig. H for b = 100, A = 1, and TV = 5000. There, we 
observe excellent agreement between the expressions for 
(5i?pRT and SE and the corresponding numerics. To verify 
the validity of Eq. (jSJ we compare the 50% probability 
width (defined as the energy width of the central r region 
that contains 50% of the probability) of P{r) to that of 
^pr.t('') finding reasonable good agreement for x < Xp^t, 
as expected. 

Does the LDOS analysis capture all the features of the 
eigenstates? The answer turns out to be negative. In 
case of the WBRM model there is still another regime 
which is not captured by the LDOS analysis. Namely, 
for X > Xoo where Xoo = b^^'^Xc the eigenstates of the 
Hamiltonian become exponentially localized. This is the 
well known Anderson (strong) localization effect. For 
X > Xoo the Hamiltonian is essentially 7i = S as if the 
diagonal E does not exist. Do we have an analogous type 
of crossover in case of the WLRM model? The answer 
must be positive because we know fl4j that in case of the 
WLRM model, the Hamiltonian H = B has multifractal 
(rather than localized) eigenstates. 



V. MULTIFRACTALITY 

The multifractal structure of the eigenstates is com- 
monly characterized by the fractal dimension D2 which 
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FIG. 4: SEpRT and 3E as a function of x for the WLRM model 
with b = 100, A = 1, and AT = 5000. The 50% probability 
width of P{r) (A) and of Pprt(»') (dashed hne) are plotted to 
verify the validity of Eq. Q for x < x^^t- SEpu.t (dotted line) 
and SB (o) were obtained from (5-Eprt = ^ X/r ^^PpRT{r) and 
SE^ = A^^r^ P{r), respectively. The numerically obtained 
5-EpRT (dotted line) is hardly visible in the plot because it lies 
on top of SE-pBT from Eq. We also show F (thin full 

line) from Eq. For this set of parameters Xc ~ 0.57 and 
Xnit ~ 5.35. 



is associated with the scaling of the inverse participation 
ratio. Given an eigenstate of Ti. which is represented by 
a column vector we define the participation ratio as 



M2 = 



(16) 



The exponent D2 is defined via the scaUng relation 



oc 



N 



D2 



(17) 



where 772 = exp((lnA/2)). The measure A/2 constitutes 
an estimate for the typical number of non-zero eigen- 
function components of the column vector. For a local- 
ized state it equals a constant number and hence D2 = 0, 
while for an extended non-fractal state it is proportional 
to the size of the matrix N and hence D2 = 1- In gen- 
eral one finds < D2 < I- The fractal dimension D2 
manifests itself in a variety of physical circumstances. 
As examples we mention the conductance distribution 
in metals the statistical properties of the spec- 

trum jl^ . the anomalous spreading of a wave-packet, the 
spatial dispersion of the diffusion coefficient and the 
anomalous scaling of delay times ■ 



N= 128 
N = 256 
N = 512 
N = 1024 
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FIG. 5: (lnA/'2(a;)) as a function of a; for = 128, 256, 512, 
and 1024. The bandwidth is fe = lOO(o), 40(0), 10(0), 4(A), 
1(<1), 0.4(v), 0.1(C>), 0.04(+), and O.Ol(x). 



In it was shown that the fluctuations of 7^2/7^2 
for Ti = B are characterized by a universal probability 
distribution. A theoretical estimation gives 



Do 



l-(7r6)"i 6>1 
2h 6 < 1 



(18) 



We notice that D2 = D{b) gives a global fit for the fractal 
dimension, where we define 



D{h) 



1 



1 + (2.346)- 



(19) 



Note that H19|) is also in agreement with the numerical 
found 20] value D2 ~ 0.7 for b = 1. For sake of later 
analysis we have found that the associated proportional- 
ity factor in Eq. (|17|l is exp{—G{b)) where 



G{b) 



1 



1 -I- (1.236)- 



(20) 



Thus, based on the above, we may say that for x = 00 
we expect to have multifractal eigenstates. We turn now 
to discuss the more general case of finite x. We want to 
see how the multifractality emerges as we increase x form 
zero to infinity. 

In the numerics we assume without loss of generality 
that the mean level spacing is A = 1. This implies that 
for 6 > 1 the threshold for mixing of levels is Xc = tt"^/^, 
while for 6 < 1 it is a;„ = (7r6)-i/2. Our main interest 
is in the non-trivial regime x > Xc- In Fig. we plot 
lnA72(a;) for several values of 6 and N. We calculate the 
average using 25% of the eigenstates at the middle of the 
spectra from a number of realizations of B summing up 
a total of 100, 000 data values for each (6, N). 

Looking at Fig. |S1 we see that for x = we have 
lnA72 = because all the eigenstates have only one com- 
ponent. We observe that this zero value persists up to a 
point X = xo{b). In principle one can argue that a;o(6) 
should be of the order Xc and (for 6 > 1) not larger than 
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F(b) ln(x/xQ)/<ln N2(oo)> 

FIG. 6: The scaled version of Fig. El {lnAf2{x)) / (lnA/'2(oo)> 
is plotted against F{b)\n{x/xo) / {lnA2(cxD)). The parametric 
crossovers at x = a;o and x = Xoo, are indicated by dashed 
lines. 



In Fig. we demonstrate this scaling. We plot 
lnA7'2(a;)/ In A7'2(oo) as a function of the scaled vari- 
able ln(2;/a;o)/ lnA7'2(oo). We see clearly the triv- 
ial crossover at x — xq and the non-trivial crossover to 
multifractality at x — Xoo ■ 

Possibly it is more instructive to describe the behavior 
of A/2 as a function of N for a given x. The interesting 
case is to have a fixed value of x which is much larger 
than Xq. As we increase N we have a multifractal growth 
N2 oc N^'^. This goes on as long as Xoc,iN;b) remains 
smaller than x. After that J72 reaches saturation, as 
implied by Eq. H23|) . The saturation value is related to 
the "width" of the LDOS, and hence has an algebraic 
dependence on the dimensional strength (x/xq) of the 
perturbation Q . It is pleasing to note that for 5 ~ 1 we 
observe F(b) « 2/3 which is related to considerations as 




^^prt ~ Vbxc- However, from Eq. H19|l it is clear that for 
b > 100 we already have D2 ~ D{oo) = 1. Therefore 
the b dependence of Xq can be neglected, and in prac- 
tice cannot be detected. We shall see below that for any 
practical purpose one can take xq « 0.15. 

As X is increased beyond xq the participation ratio M2 
becomes larger. As long as x is not too large the E term 
in the Hamiltonian dominates, and therefore the size of 
the matrix is of no importance. Indeed, we see that the 
curves in Fig. |5l are N independent for small x values. 
From this plot we find that the slope of the curves is 
given by 

-^(^^ {0.57 + 0.2 In(fe) 6>0.1 ■ ^^^^ 

For large enough x, the value of lnA72(a;) saturates to 
the X = 00 multifractal result. We call the crossover 
point Xoo- Using the knowledge of both the x < Xoo 
behavior (as described in the previous paragraph), and 
the X > Xoo behavior (saturation), we deduce that 

x^(N; b) = xoe-''^'')/PWNDib)/Fib) (32) 

Putting together all the above findings we end up with 
the following global scaling relation for the participation 
ratio: 

1 X < Xq 

(x/xo)^(^) XQ<x<x^{N-b) . (23) 

Q-G{b)^D(b) x>Xoo{N;b) 



VI. CONCLUSIONS 

We have analyzed the parametric evolution of eigen- 
states for the WLRM model. Both, the standard WBRM 
and the WLRM models, exhibit a crossover from a per- 
turbative regime where the LDOS is Lorentzian-like to a 
non-perturbative regime where the LDOS is semicircle- 
like. However, there is in both cases an additional 
crossover which is not captured by the conventional 
LDOS analysis: In the case of the standard WBRM 
model it is the well studied crossover to an Anderson 
localization regime, where the eigenstates become expo- 
nentially localized; In the case of the WLRM model it is 
the emergence of multifractality. 

We have also shown that the WLRM model possess 
the X invariance property, absent in the standard WBRM 
model and in other models that do not have a Lorenzian 
band profile. Both the non-perturbative crossover and 
the X invariance property characterize realistic quantized 
Hamiltonians. 
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